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ABSTRACT 

We study resonant CP-violating Einstein-Podolsky-Rosen correlations that may take 
place in the production and decay of unstable scalar particles at high-energy colliders. 
We show that as a consequence of unitarity and CPT invariance of the S-matrix, in 2 — > 2 
scatterings mediated by mixed scalar particles, at least three linearly independent decay 
matrices associated with the unstable scalar states are needed to obtain non-zero CP-odd 
observables that are also odd under C-conjugation. Instead, for the correlated production 
and decay of two unstable particle systems in 2 — > 4 processes, we find that only two 
independent decay matrices are sufficient to induce a net non-vanishing CP-violating 
phenomenon. As an application of this theorem, we present numerical estimates of CP 
asymmetries for the correlated production and decay of supersymmetric scalar top-anti- 
top pairs at the LHC, and demonstrate that these could reach values of order one. As a 
byproduct of our analysis, we develop a novel spinorial trace technique, which enables us 
to efficiently evaluate lengthy expressions of squared amplitudes describing the resonant 
scalar transitions. 
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1 Introduction 



In 1935, Einstein, Podolsky and Rosen (EPR) pQ proposed a remarkable experiment, 
by which the non-local nature of quantum mechanics could be tested in unstable systems 
decaying into two entangled states. Contrary to authors' initial expectations, subse- 
quent EPR-type experiments based on polarization of photons have refuted local real- 
ism, vindicating the non-local and non-causal interpretation of the quantum-mechanical 
wavefunction. The EPR paradigm has found numerous applications in modern quantum 
theory [2], such as quantum information, quantum dots and particle physics, including 
tests of CPT violation and quantum decoherence in (^-factories as predicted in certain 
theories of quantum gravity [3]. 

In this paper, we study the physical consequences of EPR entanglement on resonant 
CP violation j3] in the correlated production and decay of unstable particles at high-energy 
colliders, such as the TEVATRON and the CERN pp Large Hadron Collider (LHC) (for 
related considerations at future colliders, see jS])- As a prototype example, we consider 
the pair production of SU(3)-coloured scalar particles Si j2 via virtual gluons at the LHC, 

followed by their decay into fermions, e.g. pp — > Si^S^ — > (f f) (/'/'). For instance, 
in supersymmetric theories, such coloured scalar particles could be scalar top or bottom 
quarks, i.e. Si t2 G {t\,2 , ^1,2}, which can decay into the gauginos or Higgsinos /, /' and 
Standard Model (SM) fermions /, /'. As a consequence of cross-correlation between 
unitarity and CPT invariance, we find that the two decaying 'arms' of Si t 2 and S^ 2 are 
not independent of each other. As we will show, at least two linearly independent decay 
channels are needed to obtain non-zero CP asymmetries. The proof of this statement gets 
facilitated with the use of a new trace technique, which we introduce in order to efficiently 
evaluate lengthy scalar amplitude expressions. 

The paper is organized as follows. After this brief introduction, in Section [2] we 
review the formalism for describing 2 — )• 2 resonant scatterings, mediated by a system 
of mixed, unstable scalar particles, and present the action of CP and T transformations 
on the squared amplitudes. In Section [3] we discuss the unitarity and CPT constraints 
on CP-violating observables. Section [4] shows the connection between 2 — 2 and 2 — > 
4 scatterings, and presents numerical estimates of CP asymmetries that may occur in 
the correlated production and decay of supersymmetric scalar top-anti-top pairs. In 
particular, the scattering processes are identified for optimally testing CP violation in 
the scalar top sector of the Minimal Supersymmetric extension of the Standard Model 
(MSSM) at the LHC. Section [5] contains our conclusions. In Appendix^ we give technical 
details of a novel spinorial trace technique which we develop to analytically evaluate 
lengthy expressions of amplitudes squared describing transitions mediated by unstable 
scalar particles. 

2 Formalism for Scalar Particle Mixing 

In this section, we will first review and then extend the former formalism for unstable 
particle mixing, presented in jl]. Our aim is to study CP asymmetries in the production 
and decay of scalar particles. Motivated by supersymmetry, we consider as a prototype 
example coloured scalar particles S a , such as scalar top (£1,2) or bottom (61,2) quarks. 
We assume that these particles are sufficiently heavy, so as to decay into two fermions, 



2 




Figure 1: Scalar transitions S@ — > S a at the one-loop level. A cut on the graph according 
to the Cutkosky rules 07) / yields the absorptive contributions. 

e.g. S a -> fifj, where fa € {t, b} and f-. £ Xi,2,3,4> Xi^}- The generic Lagrangian 
describing the relevant scalar- fermion-gaugino/Higgsino interactions S a -fi-fj reads: 

^ X (cfP L + C«fP R ) f) + H.c. , (2.1) 

with Pl,r — (1 T 7s)/2. Note that explicit couplings for S a = i a ,b a in the CP- violating 
MSSM may be found, for example, in P|[TU]. 

The dynamics of the mixed S1S2 system may be described by the (unrenormalized) 
inverse scalar propagator matrix [I] 

/ s - Ml + n u (s) n 12 ( s ) \ 

An (s) = • (2.2) 

^ n 21 ( s ) s - m| + n 22 ( s ) J 

The propagator matrix arises from summing up a geometric series of the S a S/3 self-energies 
n a /3(s), as shown in Fig. [TJ The self-energies may be linearly decomposed in their disper- 
sive and absorptive parts as follows: 

= n3 P (s) + iU^(s). (2.3) 

The dispersive part of the self-energy, Il^ p (s), is UV infinite and requires renormalization. 
To this end, we consider the on-shell (OS) scheme, in which the dispersive parts satisfy 
the renormalization conditions jl] 

n d jW) = Tl d J(M*) = , jim = . (2.4) 

The advantage of the OS scheme is that around the resonant region s ~ M\ « Mf , the OS- 
renormalized, UV-finite dispersive self-energies are negligible, and only the absorptive self- 
energies Il^ s (s) become relevant. Hence, one may obtain a Born- improved approximation 
for the propagator matrix given by 



s - Mf + iUl\ s (s) ill 



abs 1 
12 



A _1 (s) = i ' • (2-5) 
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Figure 2: Transition amplitude fkfj — > fmfh via the scalars S a ,p (left), and its CP- 
conjugate amplitude (right). The grey circles indicate one-loop resummed contributions. 



Employing the standard Cutkosky cutting rules [TT] as shown in Fig. [fl we can calcu- 
late the absorptive part of the self-energies, n^ s (s). For our generic theory with scalar 

interactions S a -fi-h , we obtain at the one-loop level, 



0) = 77— A 1/2 0,m-,m?) 0(yfs - nii - m-) 



X 



167TS 

i] 13 i-J *J / 1 V % 3 l 3 %J 13 J J 



(2.6) 



where A(x, y, z) = (x — y — z) 2 — Ayz and summation over all intermediate OS states fi fj 
is implicitly assumed. Observe that the absorptive self-energies Il^ s (s) are related to the 
anti-Hermitian part of the inverse propagator matrix (12. 5p . through 



A- 1 ^) - (A~ 1 ( s )) t 
The latter is also related through the optical theorem to 



ad 



A~\s) - - (A~ 1 ( s )) t 



a 3 



X 



/ dPSx V? (^x) f 



(2.7) 



(2-* 



where the sum is over all OS intermediate states X. In addition, are the vertex 

amplitudes for the decay processes S a ^ — > X, and PSx is the respective Lorentz-invariant 
phase space. Further constraints from unitarity and CPT invariance will be discussed in 
the next section. 

It is now instructive to consider 2 — > 2 scatterings, mediated by the coloured scalar 
particles S a in the s-channel, as shown in Fig. [2j Specifically, the resonant s-dependent 
part of the amplitude A — >■ B is given by 



71 



s) = £ TO* Ao^oo vf = v3 A( S ) y A 

ct,/3 



(2.9) 



where we use matrix notation in the space of the unstable scalar particles S a in the last 
step of the equation. Moreover, the vertex functions, for the in-state A = fk fj and the 
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out-state B = f m fa, have the analytic forms 

Vf = Vl = v(PiAi)(c^*PL + C^P R )u(p k ,X k ), (2.10) 

(YbY = Krd* = u(p m Am)(C^P L + C^P R )v( Pii ,Xn), (2.11) 

where p k and X k symbolize the 4- momentum and the helicity for a given fermion f k . It now 
proves useful to define the Hermitian matrix 

Tfts) - £vj(vzr = (c^/ + c^crf)(s- m i-^) 

-2(C^/ + C^crf) mkmi , (2.12) 

for the production vertex, including an analogous definition Tj 7 for the decay vertex. 
In terms of these matrices, the squared amplitude for the s-channel resonant process 
A — > S a — > B takes on the simple form: 



\ T A^B {s)] 2 = T |a (s) Aafj{s) T ^ {s) A f a(s) = Tr 



T B {s)A{s)T A {s)A\s) .(2.13) 



Notice that the squared amplitude can be compactly written in the form of a trace over 
the scalar-particle space S a . 

For the CP-conjugate process A — > S* — > B, the transition amplitude may be 
written down as follows: 

r*-+ s = {T B y Ai P { s )v p A = v ] b a\ s )v a , (2.i4) 

where the CP-conjugate in- and out-state vertex functions are given by 

(VlT = v{p^\ h )(c^P L + C^P R )u{ Pm ,\ m ) = V§(p), (2.15) 

V A = u( Pk A k )(clfP L + ClfP R )v( Ph \i) = [Vf(p)]*, (2-16) 

with p = (E,—p). The RHSs of the last two equations give the relations of the CP- 
conjugate vertices V AB to V AB . Taking these relationships into account, we can express 

the squared amplitude of the CP-conjugate process entirely in the S^-space as follows: 
\T s A ^(s)\ 2 = £ V B dJ\7 A v\{dJ)W B (2.17) 



Tr 



T T B {s)A J {s)Tl{s){K v {s))A = Tt\t a (s)A(s)T b (s)A^s 



In fact, the last two trace expressions correspond effectively to the discrete transformations 
CP and T-reversal, respectively. Explicitly, under CP-conjugation, the resummed S a - 
propagator matrix and the production and decay matrices transform as 

CP: A(s) A(s) = A T (s) ; T A(B) (s) T A(B) (s) = Tj {B) (s) . (2.18) 

Under T-reversal, we have 

T: A(s) A\s) = A(s) ; T AB (s) T A{B) (s) = T B{A) (s) . (2.19) 

Notice that the squared amplitudes (12.131) and (12.171) are manifestly invariant under CPT 
t r ansformat ions . 
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3 Unitarity and CPT Constraints on CP Observables 



The unitarity and CPT invariance of the S-matrix introduce non-trivial constraints on CP 
asymmetries, see e.g. [I]. In particular, the equality of the s-channel forward amplitude 
A —7- A with its CPT-conjugate one A — >■ A, 



T A ^ A (s) = Tr 



s 



Ta(s) A(s) 



T s A ^ A {s) , (3.1) 



gives rise to relations among the different resonant channels in processes, such as A — > X 
and its CP-conjugate process A — > X, where X (X) is an accessible final state. To make 
such relations explicit, let us first consider the reactions of A into all possible X states. 
The sum of the squared amplitudes integrated over the phase space of the X states may 
be calculated as 

[ dPSx \T A ^ X \ 2 = J2 FSx Ti [ Ta AT * At ] = -* Tr [ T A A (A" 1 - (A" 1 )*) A f 
x J x 



= *Tr[r A (A - A^] = -2Im(7; A ^ A ), (3.2) 

where we used the relation (12.81) derived from the optical theorem. As a consequence of 
the CPT constraint (j3.1j) . the total unpolarised cross section of A to all possible X is 
equal to the corresponding unpolarised cross section of the C-conjugate states A to all 
possible C-conjugate states X. 

Assuming that the dominant source of CP violation is due to scalar-particle mixing, 
we define CP-violating observables pertinent to a resonant reaction A — > S a — > B and its 
CP-conjugate one A ->■ S* -»■ B: 

A™ B {s) = J (\T A ^ B (s)\ 2 -\T A ^ B (s)\ 2 )dPS B . (3.3) 

Since we sum over all particle helicities and integrate over the P-invariant phase space, 
the CP-odd observable A^ B is also odd under C, but even under P transformations. 

We will now show that thanks to unitarity and CPT invariance, at least three 
linearly independent decay channels of S a are required in 2 — > 2 scatterings, in order to 
obtain non-zero C/CP-odd observables, such as A^ p >B (s) given in (13. 3p . Our proof for 
a two-particle scalar mixing system will then be generalized to mixing systems with n 
scalar particles. In Section m we extend this theorem to the correlated production and 
decay of two unstable-particle systems in 2 — > 4 scatterings, where we show that only two 
independent decay channels are sufficient to obtain a non-zero C and CP asymmetry. 



3.1 CP Conditions in Two-Particle-Mixing Transitions 

We will now show that the C- and CP-odd observables A < ^ B (s) defined in (13. 3p vanish, 
if only two independent decay channels are open in 2 — > 2 resonant scatterings A — > S a — > 
X. To this end, let us assume that X = A', B, where A' is the same two-particle state as 
the initial state A but has different momentum configurations. As a consequence of the 
CPT-invariance and unitarity relations in (13.11) and (13.21) . we have the equality 



\r A ^ A '\ 2 + I dPs B \r A ^ B \ 2 = I dPs A ,\T A ^ A '\ 2 + I dPs B \r s A ^ B \ 2 . 

(3.4) 
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In terms of CP observables, the above equality can be rewritten as 

= J dPS B (\T S A ^ B \ 2 - \T S A ^\ 2 ] = - [ dPS A , ( \r^ A '\ 2 - \T^ A '\ 2 



= -A^ A (a). (3.5) 

Since the decay channel of S a into states A' is not kinematically independent to the 
production channel A — Y S a (they just differ by an overall phase-space factor), we have 
that = 0, as a consequence of CPT invariance, implying also that A A ^ B (s) = 0. 

Hence, we can conclude that more than two linearly independent decay channels would 
be needed to obtain non-zero C and CP asymmetries in 2 — > 2 resonant scatterings. 

The above conclusion can also be obtained by an explicit calculation of the squared 
amplitude A — > S a — > B, in a two-particle mixing system; the generalization to a system 
with n particles follows in the next subsection. Our calculation gets facilitated by a 
spinorial trace technique, where the 2x2 production and decay matrices T a ,b an d the 
inverse propagator matrix A" 1 are expanded in terms of the two-by-two unit matrix 
1 2 = cr° and the three Pauli matrices a 1,2 ' 3 . Technical details of the trace technique are 
given in Appendix [A] In detail, we obtain for the squared amplitude 

rnH A —1,1/ rpp I A -1,A\* 

\r^ B \ 2 = Tr \T A A T b A*] = | (A 1 )| 2 Tr ^^ > ( 3 - 6 ) 

where the expansion coefficients are given in Appendix IA.2I Since the CP-conjugate 
amplitude is obtained by the interchange of the Hermitian vertex matrices T A ^ T B , the 
only contribution to the CP observable A A ^, B (s) comes from the fi <-> p antisymmetric 
parts of the amplitude squared. These are the terms proportional to the the Levi-Civita 
tensor e^px contained in the trace expression over the generalized Pauli spinors on the 
RHS of (13. 6p [cf. (lA.lOp ], Specifically, up to overall phase-space and other kinematical 
factors, the CP observable A A ^ B (s) is 

A^p B ( s ) oc e, upX T% T B Im( A _1,p ) Re(A~ 1 ' A ) . (3.7) 

The 4-vector Im(A~ 1,p ) derives from absorptive parts of the self-energy, 

Im(A-^) = l - Tr [^U ahs ] = 1 (pS a T^ + PS B T^ , (3.8) 

where PS^^ denote phase-space factors associated with the A and B states. Evidently, 
three linearly independent channels or 4- vectors in the generalized Pauli-spinor space are 
needed, e.g. T^ BC , in order to obtain a non-zero CP asymmetry A A ^ B in (13.71) . 

3.2 Beyond the Two-Particle Mixing 

The above result obtained for a two-particle-mixing system can now be generalized to 
mixing systems with n unstable particles S a , with a = 1, 2, . . . , n. 

Our starting point is the n x n propagator matrix A given by 

A(s) = [D" 1 ^) + ill(s)]- 1 . (3.9) 

The n x n propagator matrix A(s) consists of a Hermitian matrix D~^(s) = 5 a/ 3(s — 
M 2 ) + n^ p (s), which may include dispersive contributions, and a general anti-Hermitian 
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matrix z'II Q| g(s) = iU^p(s) that describes the absorptive effects. Dropping the explicit 
s-dependence, it can be rewritten in the more convenient form: 

A = (D — z D n D) [l + (n D) 2 ] _1 = [l + (D n) 2 ] 1 (D — i D n D) , (3.10) 

where 1 is the n-dimensional unit matrix. 

The cross sections for the 2 — > 2 resonant process A — > S a — > B and its CP-conjugate 
one A — > S* — > B may straightforwardly be evaluated by the trace expressions 



a = Tr [T B A T A A f ] , a CP = Tr [T A A T B A f ] 



(3.H) 



where phase-space integration is implied. From these, the CP-violating difference of the 
two cross sections, A<t cp = a — <x CP , can be defined, which is calculated as 



CP 



- Tr 

2 



Tb (A — A*) T A (A + At) - T B (A + At)T A (A-At) 



= -2i Tr | Tb D n [1 + (D n) 2 ]" 1 BT A [t + (D n) 2 ]" 1 D 

-T B [1 + (Dn) 2 ]" 1 DT A Dn[l + (DlTj^Dj . (3.12) 

From the last expression in f l3.12p . we observe that the CP- violating difference Aa CP can 
be cast into the compact form: 

Aa CP = -2zTr{DII (l + (D n) 2 )" 1 DT A , (l + (Dn) 2 )"' DT B }, (3.13) 

where [X , Y] = XY — YX is the commutator for two n x n matrices X and Y . Clearly, 
if Tb = cT A , with c G C, the commutator in C I3 . 1 3[) vanishes. Further, from the cyclicity 
of the trace, 



Tr{A[B,C]} = Tr{C[A,B}} = Tr{B[C,A}}, 
we have that the CP-violating difference vanishes, i.e. A<r CP = 0, if 



or if 



(l + (Dn) 2 ) 1 D T A , Dn 



(l + (Dn) 2 ) 1 D Tb , Dn 







. 



(3.14) 



(3.15) 



(3.16) 



Since the absorptive contributions can in general be parametrized by the different open 
decay channels, n = aT A + (3Tb + 7Tc, then Acx CP will vanish if Tq = (only 2 open 
decay channels), or if Tc = c\T a + c{Tb is linearly dependent. Consequently, at least three 
linearly independent decay channels are needed to construct non- vanishing CP asymme- 
tries in general 2 — > 2 scattering processes. As we will see in the next section, the validity 
of this theorem may be extended to unstable states which are pair produced and decay in 
correlation in 2 — > 4 scattering reactions, where however only two independent channels 
are needed. 

We conclude this section by remarking that in the absence of a particle mixing, 
the decay matrices T AjB ,c an d n are diagonal, and the commutator in (I3.13P vanishes, 



S 




Figure 3: Correlated production and decay of S a particles in 2 — y 4 scatterings at high- 
energy hadron colliders. 

leading to Act cp = 0. This corresponds to the case where each unstable particle couples to 
a disjoint set of states. This situation can happen naturally, only if the unstable particles 
carry their own conserved charges. In minimal supersymmetric theories, the mixing of 
light squarks or sleptons is suppressed by their masses, such that CP asymmetries are 
sizable, only if particular resonant conditions are met [UEI21EE2] • As we wm see i n t ne next 
section, we find sizable CP asymmetries in the strongly mixed stop or sbottom sectors for 
a wide range of MSSM parameters. 



4 EPR Correlated Production and Decay of Scalars 



In this section we study the EPR-correlated production and decay of unstable scalar 
particles at high-energy colliders. Specifically, we consider the resonant part of 2 — y 4 
processes, mediated by coloured or charged scalar particles S a , as displayed in Fig. [3j 

The transition amplitude for the 2 — y 4 partonic process qq —y S a S^ — y (/ /) (/' /') 
has the analytical form 

r^\h,k 2 ) = J 9 (k 1 ,k 2 )Vl(k 1 )A(s 1 )A(s 2 )V B (k 2 ), (4.1) 

where the initial gluon current factorizes in the amplitude, J g (ki, k 2 ) = {k\ — k 2 )^ Q 
Absorbing all kinematical factors into an overall normalization constant A/", the cross 
section for the 2 — y 4 scattering is given by 



AfTi 



T A ( Sl ) A( Sl ) A(s 2 ) T B (s 2 ) A\s 2 ) A\ Sl ) 



(4.2) 



where s = (k\ + k 2 ) 2 3> M 2 , si t2 = k\ 2 and phase-space integration is understood. 

Since we are interested in the dominant part of the cross section, we expand cr 2 ^ 4 (s) 
about the resonant region y/^ = (Mi + M 2 )/2. Hence, we obtain for the resonant (res) 
part of the cross section 



AfTi 



T A (s) (A(s)f T B (s) (At( s -)) ; 



(4.3) 



3 Here we discuss the LHC suppressed, but more intuitive one-gluon exchange process. Other processes, 
such as gg — > S a S^ — y (//) (/' /'), exhibit similar analytical features and can be treated in analogous 
fashion. 
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where we have ignored sub-dominant terms proportional to (s\ — s-i)js. 

To linearize the expression H4.3[) in terms of the propagator matrix A(s) and its 
Hermitian conjugate A^(s), we define 



A(x;s) = [xt-M' + ilL^is)]- 1 
with A(s; s) = A(s), such that 

(A( S )) 2 = (A(,;s)) 2 = 



^ A / \ 



(4.4) 



(4.5) 



With the help of this last relation, the resonant part of the 2 — > 4 cross section can now 
be recast into the form: 



a 



2->4 



N ° 2 T 
<9x cfy 



T A (s)A(x; s)T fl (s)A%;s) 



(4.6) 



x=y=s 



Note that this form is very similar to the squared amplitude of a 2 — > 2 resonant process 
[cf. (I3.6|) ]. For the specific two particle mixing case, the resonant part of the 2 — y 4 cross 
section is given by 



a 



2->4 



TV 



O 2 



-1,1/ 1 



.r 



9x9?/ detfA" 1 ^; s)] detfA" 1 ^; s)] H 



x T^Tg Tr [a^a u a p ax} . (4.7) 



From this last expression, we may calculate the CP-odd part of the 2 — y 4 cross section, 
which may conveniently be written down as 



2->4\QP 



AMe, upX T»T» B lm 



d A~ 1,/3 ( 



x; sj 



<9x det[A _1 (x; s)] 



Re 



9 A~ 1,A i 



x; s 



dx det[A _1 (x; s)] 



(4- 



where x is set to s after differentiation. Unlike in the 2 -> 2 scatterings [cf. (13.71) ]. we 
may convince ourselves that for 2 — y 4 scatterings, only two linearly independent decay 
channels would be sufficient, i.e. 



Im (A- 1 -'^)) = aT p A + /3T 



B ' 



(4.9) 



in order to obtain a non-zero CP-odd contribution to the cross section. Specifically, it is 
not difficult to see that the CP-odd part (cr^ i ) < ^ P contains terms proportional to 



rp/I rpy 

tfivpX 1 A 1 B 



d Re A- 1 '' I 



x; s 



dx 



ReA" 1 ' A (s) Imfdet [A(s)] 



(4.10) 



which are clearly non-zero, if only two channels as given in (I4.9j) are assumed. 

In our discussion, we have ignored other one-loop effects, such as gluon-scalar-scalar 
vertex corrections. These corrections are sub-dominant, but become relevant in restoring 
the gauge-fixing-parameter independence of the gluon propagator by means of Ward- 
Takahashi identities (for a recent review, see [33]). In particular, in the Feynman gauge, 
one can show that all these effects are non-resonant and can thus be safely neglected. 
A detailed study of these higher-order effects is beyond the scope of the present paper. 
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4.1 Application to Stop CP Violation at the LHC 



As an application of our results, we consider the pair production and decay of super- 
symmetric scalar top (or bottom) quarks [5J. These coloured particles can be copiously 
produced via partonic QCD interactions at hadron colliders. Their sizable mixing can 
lead to large asymmetries, due to non-vanishing SUSY CP phases in the heavy squark 
sector q, q = t,b. Their dominant decay modes are 



?1,2 
?1,2 
91,2 



— >■ 



Q + 9 , 

1 + Xl,2,3,4 > 



<I+ Xl,2 ■ 



(4.11) 
(4.12) 
(4.13) 



Instead, the decays of q\p into neutral bosons are closed for sufficiently small q\-q 2 mass 
splittings, i.e. 



q 2 t4- gi + iZj 2)3 , 
q 2 q! + Z°. 



(4.14) 
(4.15) 



Likewise, their decays into charged bosons are closed for sufficiently small stop and sbot- 
tom mass differences, i.e. 



?1,2 ~fr Ql,2 + H 
91,2 7* q'l,2 + W ~ 



(4.16) 
(4.17) 



Stop (and sbottom) decays at the one-loop level within the CP-violating MSSM have been 
calculated in [7], and in the CP-conserving MSSM in [8]. 



For our numerical example, we choose the MSSM parameters Mq = Mq = Mi 



470 GeV, Mr 



D — 

M£ = 120 GeV, \A t \ = \A b \ = \A T \ = 110 GeV, M 2 = 2M l = \fi\ = 
250 GeV, nig = 710 GeV, tan/3 = 5, and the phases <pA t — 7r /4, 4>A b = 4>a t = <j>Mi — 
4>fj, — 'Pg — 0. At tree level, the resulting stop masses, widths, and mixing angle are 

mz 



^{1,2} 

are m-± 

X {1,2} 



{485,514} GeV, T i{i2} = {5.3,9.1} GeV, cos^~ = -0.725. The gaugino masses 
{187,321} GeV, m^o i 2 3 4} = {115,194,256,322} GeV. 

As an illustrative example, we study stop pair production via the process qq — > 
i a i*, decaying into charginos: t a — > bxf (A), £* — > bx^ (B). At the resonant region, 
(m^ + mi 2 )/2, the transition amplitude squared and its CP-conjugate one follow 



from (I4.3p to be 



\T(s)\ 2 

\m\ 2 



Tr 
Tr 



(4.18) 
(4.19) 



T A (s) (A(s)) 2 T B (s) (At(s)) 

T B (s) (A(s)) 2 T A (s) (At(g)) 
In Fig. HI we show the phase <pA t an d the stop mass dependence of the CP asymmetry 

m 2 -m 2 



A 



|7T+|7T 



(4.20) 



We observe that the CP asymmetry can reach values of order one, thanks to a maximal 
mixing of the two stop states in production and decay. The CP asymmetry would decrease 
as stop mass splitting increases, which can occur for large mass differences Mq — Mfj. 
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Figure 4: Tae CP asymmetry A, see ^4.2(J\ ), for stop pair production qq — » £ a (s) i* (s) and 
decay t Q &Xi" ; t* a 0X2", at V^" = {m^ + rrif 2 )/2, as functions of 4>A t an d m^. Also 
shown are the branching ratios BR(ti bxt) (solid, blue), BR(ti — > bxt ) (dashed, red), 
BR(i 2 — > bxt) (dotted, black), BR(t 2 —> bxt) (dash-dotted, green), as functions of 4>A t 
andm^. In the m^-plots (c) and (d), Mq = Mjj = is varied, whilst the other MSSM 



parameters are fixed as given in the text of sub section U7I 



In Fig. [51 we study in more detail the dependence of the amplitudes [cf. (14.2 



|T( S i,s 2 )| 2 

ir( Sl ,s 2 )| 2 



Tr 
Tr 



T A ( Sl ) A(si) A(s 2 ) T B (s 2 ) At( S2 ) A+(si) 
T B (si) A( Sl ) A(a 2 ) T A (s 2 ) ti(s 2 ) tf( 8l ) 



and the corresponding CP asymmetry as defined in (I4.20p . as functions of 5^ 
/ s 2 )/2. In our numerical estimates, we use the parameterization 



(4.21) 
(4.22) 



(4.23) 
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Figure 5: The dependence on = (y/si — y/s2)/2 of (a) the amplitude squared 

|7~(si, S2) | 2 \4-%l\) (up to an arbitrary normalization Af), with ^/si,2 = V^" ± 5y/s 

mf )/2, for stop pair production qq — > t a (si) t„(s 2 ) and 



mi 



the resonance region 

decay t a — > bxt, t* a ~^bX2 (solid, blue), and the CP-conjugate amplitude squared \4-2£fy 
for the process t a —> bxt > ~^bXx (dashed, red), and (b) the corresponding CP asym- 
metry ^TJjUj) , with MSSM parameters as given in the text of subsection \4-l\ 



in* 

"AT" 0.5 




(a) 



-15 15 

S^J [GeV] 




(b) 



-15 15 
5^/1 [GeV] 



30 



Figure 6: The same as in Fig. [5] /or (a) the stop pair production qq — > t* (#2) an d 

decay t a — > bxt, ta^bXi (solid, blue), and the CP-conjugate amplitude squared for the 
process i a — > bxt > ^0.^^X1 (dashed, red), and for (b) the corresponding CP asymmetry. 



around the resonance region yfs = (m^ + m t ~ 2 )/2. As can be seen from Fig. |5l the CP 
asymmetry A can be of order 1 for the experimentally testable interval of 5^ ~ 20 GeV. 

Finally, in Fig. E]we display numerical estimates of stop decays into the same states: 
t a bxt C^)) i* a — > t>Xi (^)- The corresponding CP asymmetry A is non- vanishing, 
of order 5%, only for non-zero values of 5y/s. Hence, CP violation in the same channel 
would be more difficult to probe for the stop scenario of the MSSM under study. 
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Before closing this section, it is interesting to comment on the discovery potential 
of the LHC in measuring a non-zero CP asymmetry. As shown in [15], the minimal 
luminosity £ to observe a non-vanishing CP asymmetry A may be estimated by 



where n is the number of standard deviations (confidence level) and a is the cross section 
of the signal events. Obviously, to obtain realistic values of statistical significances and 
luminosities requires detailed Monte-Carlo simulations. Since such a study is beyond the 
scope of the present paper, here we only offer an initial estimate. 

For our benchmark scenario, with mi « 500 GeV, the cross section for stop pair 
production is cr p (pp — > t a t a ) ~ 0(700) fb, for a = 1,2, at the next-to-leading order, 
according to the public code PR0SPIN0 [16J. The stop branching ratios into charginos are 
found approximately to be BR(t a — > bxf) ~ 20%, for each j — 1,2 and a = 1,2 (see 
panels (b) and (d) in Fig. H]). Correspondingly, the leptonic (£ = e,/x) branching ratios of 
the charginos are BR(xf -+v t f)™ 30%, using the formulae quoted in [17] and the fact 
that the sneutrinos i>n decay invisibly in our scenario, i.e. U£ — > v^Xi- 

Taking all the above estimates into account, we find a net signal cross section a ~ 
0(3) fb, for stop pair production and decay. In addition to using standard cuts to reduce 
the background and isolate the signal, one has still to identify the charges of the final 
leptons and the mother chargino states, Xi or Xii f° r a proper measurement of the 
CP asymmetries. Recently, a similar study of triple-product asymmetries in stop pair- 
production and decay [T5] has shown that the method of kinematic event selection allows 
such an identification of the decay chains, if the particle masses are known, presumably 
from other measurements. Because of the kinematic cuts and the event selection criteria, 
the signal events are reduced typically by an order of magnitude, with the surviving 
CP-even background being similar in size to the signal [18]. Based on these results, our 
CP-asymmetry is expected to be reduced by half, such that we are left with A ~ 50% in 
the best case scenario. Inserting these numbers into Eq. (I4.24p . we find C ~ n 2 10 fb -1 . 
Thus, LHC luminosities of order a few 10 fb~ would be sufficient to observe resonant 
EPR CP asymmetries of order one. 




(4.24) 
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5 Conclusions 



We have analyzed resonant CP-violating Einstein-Podolsky-Rosen (EPR) correlations 
in the pair production and decay of unstable scalar particles at high-energy colliders. 
We have shown that as a consequence of unitarity and CPT invariance of the S-matrix, 
a minimum number of linearly independent decay matrices associated with the unstable 
scalar states are necessary in order to obtain a non- vanishing result for CP-odd observables 
that are also odd under C-conjugation, but even under P transformations. For 2—7-2 
scatterings, at least three such independent decay matrices are needed, whereas for 2—7-4 
scatterings, only two independent decay matrices are required. Even though our study 
involved the mixing of scalar particles only, we may safely conjecture that our results will 
also hold true for the mixing of unstable fermions and vector bosons. 

As a direct application of the above theorem, we have presented numerical estimates 
of CP asymmetries for the correlated production and decay of supersymmetric scalar top- 
anti-top pairs at the LHC. We have explicitly demonstrated that the CP asymmetries 
could be sizeable, reaching values of order one and so making them directly testable. 
Our analysis also led us to the development a novel spinorial trace formalism, which 
helped us to efficiently evaluate lengthy expressions of squared amplitudes describing 
resonant scalar transitions. The formalism presented here could be extended to other 
unstable-particle systems, predicted in minimal extensions of the Standard Model that 
include mixing of Higgs-boson states [19,20J or the mixing of heavy unstable neutrinos [2Tj . 
It is interesting to study systematically possible phenomena of resonant CP-violating EPR 
correlations in these systems both for the LHC and future high-energy colliders. 
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A Spinorial Trace Technique for Two Scalar Mixing 

Here we present our conventions for our spinorial trace technique, as well as useful iden- 
tities and analytical expressions for the case of the two-unstable-particle mixing. 

A.l Preliminaries 

A basic element of our formalism are the generalized Pauli matrices [6]: 

°» = (oV), ^ = (o ,-*), (A.l) 

with 



x o 1 y ' V 10 / V * / V _1 , 

Note that these are also the generators of the SL(2, C) group. Then, any 2 x 2-dimensional 
complex matrix A can be expanded as 

A = AX = A»a» , A- 1 = ^ , (A.3) 

where A V A V = detA and the complex components A u may be evaluated as 

A v = -Tr[Aa v ] . (A.4) 

In calculating traces, the following relations are useful [6]: 

a"a u a p + a p a v a" = 2(g^a p - g^a v + , (A.5) 

/aV - /aV = -2ie IH ' pX a x , (A.6) 

aVff p + ff p aV = 2(^^ - g^d" + g'Pa") , (A.7) 

- ff p (7"ff" = 2i£" vpx a x , (A.8) 

with the convention £0123 = +1- For instance, the traces involving two and four generalized 
Pauli matrices are given by 

Tr [a"a v ] = 2g pv , (A.9) 
Ti[a p a v a p a x ] = 2 (g^g px - g pp g vX + g vp g px - ie pvpx ) . (A.10) 

A.2 Expansion of the Scalar Propagator Matrix 

We can now use the 4-dimensional basis of the generalized Pauli matrices cx M to expand 
the inverse propagator in ()2.5|) . i.e. 

s - M? -ImII?!? s (s) \ / II?? 8 (s) Reliefs) 
A'Hs) = }+i \ . (A.ir 

Imn^ s (s) s - Mf J \ ReHfJf( s ) IB%*( N 
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To this end, we use f lA.3j) to write the 2x2 propagator matrix A(s) as 

A(s) 



A" 1 ! 



s cr^ 



detfA-^s)] ' 

where the complex 4- vector components A~ 1 (s) are given by 



(A.12) 



A^(s) 
A^(s) 



\{Ml + Ml) + \ L 



?Ren* s (s) 
■ilmn^(s) 



22 



\{Ml-Ml) + I 



(A.13) 

(A.14) 
(A.15) 

(A.16) 



and 



det [A' 1 (a)] 



{s-Mt){s-M*) + \Ii\f(s 



irff( s ) ntf( s ) 



+ z 



IT^( S ) (a - Mf) + Ut s {s) (s - Ml 



abs / 



(A.17) 



Correspondingly, the Hermitian decay matrices Ta,b,c for the final states A = (fkfi), 
see (12.121) . etc., may be each decomposed as T(s) = T M (s) cx M , with the real coefficients 



Us) 
Ti(s) 

Us) 



1 r 

2 



T n (s)+T z ' A (s 



22, 



ReT 12 (s) 
-ImT 12 (s) 



T 11 (s)-T 22 (s) 



(A.18) 

(A.19) 
(A.20) 

(A.21) 



where 



T 12 (s) 



(|CS L T + ICS""! 2 ) (» - mj - mf) - 4Re (c^f ) m t m r . 



Note that the partial scalar decay widths at tree level are 



16ttM3 1 a> 



(A.22) 



r Q ( s a — > f k fj 



(A.23) 
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